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ABSTRACT 

Path integral solutions with kinetic coupling potentials oc p\pi are evaluated. As examples I give a Morse 
oscillator, i.e., a model in molecular physics, and the double pendulum in the harmònic approximation. 
The former is solved by some well-known path integral techniques, whereas the latter by an afhne trans- 
formation. 



1 Introduction 



Quantum mechanics is about physics on the atòmic level, i.e., wave-functions and energy levels 
of àtoms and molècules. The simplest quantum mechanical system is the harmònic oscillator 
with the characteristic h(n + 5) spectrum [|]]— j|] which was also the first system solved by means 
of the path integral. The important Coulomb system, i.e., the hydrogen atom, was only solved 
later on by Duru and Kleinert || H|. In this Paper I want to demonstrate the technique of 
path integration for some particular models with kinetic coupling. The first example consists of 
the potential force of two molècules, where the anharmonic interaction force is described by the 
Morse potential with an additional kinetic coupling oc P1P2, respectively oc X1X2, the indices 1 
and 2 referring to particle one and two. The second example is the double pendulum, a system 
which is coupled by a oc <pnp2 term. In the harmònic approximation the double pendulum is 
described by two coupled oscillators. In fact, these two examples are not too difhcult to solve, 
but they demonstrate in a nice way the applicability of path integration techniques for systems 
which are coupled by means of the kinetic energy. Such problems seem not to have been studied 
by path integrals yet, whereas systems which are coupled through only a coordinate dependence, 
i.e. oc X1X2, can be found in e.g. ||. 



2 Morse Oscillator with Kinetic Coupling 



The first example is a model in molecular physics, which takes into account the anharmonicity 
of the inter-molecular potential (e.g., H2O, H2S, C2H2), by means of the Morse potential 0-0 



V(x x -x 2 ) = \( 



1 — a e" 



-(3(x 1 -x 2 -r ) 



(1) 



where xi,x 2 denote the position of the first and second particle one with masses mi,m.2, ro the 
equilibrium position, and a, /3, X are positive parameters. h denotes Planck's constant divided 
by 2tt. The Hamiltonian of this system for the two partides has the form (with the momentum 
operators P\^ = (h/ï)d xi2 ) 



H = —P? + —P? + kPiP 2 + -ae^' 11 - 12 ^ 
2mi 2m2 \ 



(2) 



The corresponding Hamiltonian path integral for the time evolution from £1,2 (í') to xi^it") in 
the time interval T = t" — t' can thus be written as [01 



x 1 {t")=x'{ 



x 2 (t")=x% 



K{x'{,x , 1 ,4,x' 2 ;T) = J V Xl (t) J VP^t) J Vx 2 (t) J ' VP 2 (t) 

xi(t')=x' 1 X2(t')=x' 2 



x exp 



Pl±l + P2X2 
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-Pt 



2m 1 1 2m 2 2 
—KP1P2 - A(l - ae-^ 1 -^ ^ 

We can switch to center-of-mass X and relative coordinates x r by means of 
X = [iixi + 112X2 , x r = xi - x 2 - r , 

P = P 1+ P 2 , P r = /í 2 Pi - /xi P 2 , 



dt 



(3) 



(4) 
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where (i = 1, 2) /íj = rrii/M, (ni + \i 2 = 1), M = mi + tti-2 is the total mass, and ji = m\m 2 jM 
is the reduced mass. Vice versa we have 

Xi = X + —X r , x 2 = X - —x r , Pl = P T + /iiP , P 2 = fi 2 P ~ Pr ■ (5) 
The Hamiltonian path integral in center-of-mass and relative coordinates is therefore given by 

X(t")=X" x r {t")=x'l 

K(X",X',x",x' r ;T) = J VX(t)JvP{t) J Vx r {t) J VP r {t) 
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Xf^t ^ — X j. 



x exp 



PX + P r x r 



(- 



\2M 



+ K/Ul/X 2 )P 2 



í(/í2 — Hi)PPr — A^l — ae 



2^1 



(ix, 
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dt 



(6) 



We can perform the momenta integrations to obtain the Lagrangian path integral for the coupled 
system 



X(t")=X'' 



X(t')=X' 



X 7' ^ t J X j. 

K(X",X',x'ï,x' r ;T)= J VX(t) J Vx r (t) 

Xyi^b ^ Xj, 
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and I have used the abbreviations 
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1 K, 

a= 2M + Kl·i2l·l2 ' d = 2~il ~ K ' b= 2^ 2 ~ ^ 



(7) 



(8) 



and D = ad — b 2 = ^(l/mim 2 — k 2 ). The limiting case k = is easily recovered. The lattice 
formulation of the path integral (0) is given in the usual way by 



K(X",X',x"x' r ;T) = lim ( —rl 
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nL dX >L dx ' 
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(9) 



Here we have taken as usual e = T/N = (t" — t') /N in the time slicing, Xj = X(tj),x rj = x r (tj) 



(j = 1,...,N), AXj = Xj -Xj-i,Ax 

r 



íCr,_i, as N — > oo,e — > 0, T fixed, and 



X" = X(t"),X' = X(t'),x" = x r (t"),x' r = x r (t'). The path integral in the center-of-mass 
coordinate is coupled to the relative coordinate x r via x r X, thus the X-path integration is that 
of a free particle with a magnètic field, which is time dependent but not X-path dependent. The 
classical solution for the center-of-mass coordinate X is found to be 



x c1 (í) = x T x (t - i') - ^ fx r {t)dt , 



with corresponding classical action 
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We can therefore separate off the center-of-mass path integration X according to [||, 11 

'i d (X"-X'f\ 
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(13) 



Again, the case k = gives complete decoupling (note 1/Ad = fi/2 for k = 0). 

The path integral in the second line of (13) is that of the Morse potential which can be found 
in the literature S. Setting 
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(14) 
(15) 



where Gk is the (energy-dependent) Green function corresponding to the kernel Kk- We find by 
means of the space-time transformation (Duru-Kleinert transformation) [|, ||] for the Green 
function 

2 T _i\ 
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Here I have abbreviated 
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and denotes the larger/smaller of x' r ,x". From the poles of Gk the bound states are given 
by 



h 2 K 2 1 - (íMk ( 1 



E,, 



2M 1-2/í/í \2(jl 



— - n)h 2 f3 2 (n + i) 2 + 2h(3{n + \)\X[ — -K 
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(18) 



The first term is just the continuous spectrum which is not relevant for the discrete spectrum.The 
bound state wave-functions are determined by the residua of the poles of Gk, and we find 

§ 0, HI 

é KX //3(2£-2n- l)n!(2a0 2Ç - 2n ~ 1 



^n íK (X, x r 
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(19) 
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The continuous spectrum is determined by the cut in Gk and I obtain 



*kAX^r) = ^ r(i* - Ç + \) e-^W^{2a^) , (20) 



with the energy-spectrum 

h 2 K 2 l-uMK / 1 \ 

E * = ^T^ + {m + ^) hk '■ < 21 > 

3 The Double Pendulum 

As the next example we consider the double pendulum with two masses mi,rri2 attached to two 
strings with lengths h,l2- g denotes the gravitational acceleration. Its Lagrangian is given by 

m 1 + m 2 2 -2 , m 2 ,2 -2 

£ = 2 1 ^ + ~2~ 2<f2 

+m2hhfii&2COs(tpi — (f2) + (mi + m2)ghcosipi + m2gl2COSip2 ■ (22) 

In the harmònic approximation we have (1% = I2 = l, M = m\ + m-2) 

M 2 - 2 1 m 2 ; 2 -2 1 ,2 ■ ■ M , 2 m 2 ; 2 / OQ A 

C = Y {f>x ~2 ^ 2 ipiíf2 ~ T 3 ^ ~~ T" 992 ' ' ^ 

which I consider in the sequel. The path integral formulation for the latter case has the form 

¥>l(f)=<pi <P2(V)=<f>2 



x exp 



M ,2 -2 m 2 ,2 -2 ,2 ■ ■ ^ , 2 m 2 , 2 

y» </>i + ~2 1 f2 + m 2l fW2 ~ y 9 l( Pl - ~2 ™ 



dt\ . (24) 



The path integral (^4|) is of the form of a path integral for a general quadratic Lagrangian. Its 
solution is given by, e.g. |, |, g, §, ||] 



*(rf,rf,rf,*4;r) = (^'Y^ ( - (>(*",*')) . (») 

where I have denoted = (yi,^), and Sç\{ip" ,ip') is the corresponding classical action of the 
system. The determinant appearing in ( |25|) is called the Morette-Van Hove determinant g, [lq| . 
In order to determine the classical action, one considers the classical system of the coupled 
differential equations and switches form the coordinates (^1,^2) to normal coordinates (£1,^2) 
by means of an affine transformation flíq ]. This is done in the following way: We write 

£ = ±£A£* - l<pK<ffi , (26) 

where A, K are n x n matrices describing the classical system with its potentials and couplings. 
The classical system is diagonalized by considering the eigenvalue equation 

det(K - lü 2 A) = , (27) 

with eigenvalues uJ 2 . This is performed by a matrix C, the affine transformation. The corre- 
sponding differential equations in the normal coordinates £ = C<p (k = 1, . . . , n) are 

Í k + = , (28) 
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with the uij. the eigenfrequencies corresponding to the kth component. In our case we have n = 2, 
and we find 



= yy (1 - r) (yi® ~ ^a(*)) , £ 2 (í) = y^(l + r) (</>i (í) + r</> 2 



(29) 



where r = y/m^/M, uip = (Mg/m\ï){l ± r). Vice versa we have 



, V2(<) 



1 



-6(í) + 



í\/2M\\/l-r VT+Ty ' rzv " y /\/2M VrVl - r ' ry/T+r 
Therefore we can write down the solution of the path integral (p4j) in the following form Q- 

= n 



(30) 
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m l w l /-, n// / /x2 , / // t/\2\ 

m2UJ2 (l + r)(K + V 2 ) 2 + K' + r^) 2 ) 



sin u^T 

-a/Ï - r 2 ((^'i - rtp 2 )((p'( - rtp'í) + (^ + rtp' 2 )((p'l + r^jf) 

II E KM)*nM>- i[ " k+x,2)T . 

fc=l,2n fc =0 



The energy spectrum has the form 



n.i,n2 



^(üJiTli + üJ 2 "-2 + 1) = h 



Mg ( 11712 , 1 

- Ini + n 2 + A / — (m - n 2 ) ] + 1 



mi 



.1/ 



and the wave-functions are given by (^ as in (pS|)) 

1/4 " _ 



*n fc (¥>fc) 



1 



2 n fcn fc ! 



exp 



m k uj k 2 \ / /m fc Wfc 



The expansion (|33| ) can be derived by means of the Hille-Hardy formula y, 17] 



(31) 



(32) 
(33) 



(34) 



(35) 



4 Conclusions 

The results clearly show the applicability of path integration techniques in systems with kinetic 
energy couplings. The first example of a Morse potential interaction, a model known in molecular 
physics, was transformed in center-of-mass X and relative coordinates x r , followed by the path 
integration of the center-of-mass coordinate, thus achieving a complete separation of variables in 
X and x r . This left the remaining path integral as a path integral Kr of the Morse potential, and 
the corresponding Green function Gk could be written down. The usual analysis of Gk gave the 
wave-functions and the energy spectrum of the discrete and continuous spectrum, respectively. 
The bound state energy spectrum has the typical feature of the energy spectrum of a Morse 



5 



potential, however modified by the coupling constant k of the KPiP2-teim. Let us note that the 
choice of the Morse potential for the interaction in the relative coordinate is but an example. We 
can take any potential V(x r ) = V{\x\ — a^D- This has the consequence that in (|ï~3|) the potential 
term Xa 2 (e~ 2 ^ Xr — (2/a)e~ l3xr ) is replaced by V(x r ), and the corresponding path integral is 
exactly solvable for any exactly solvable one-dimensional potential, modified by the mass-term 
1/M. 

The double pendulum in the harmònic approximation could be separated in terms of normal 
coordinates £ which are uniquely determined by the affine transformation C in terms of the 
original coordinates <p. The path integral can be stated explicitly, and the wave-functions and 
the energy spectrum as well. The formulation in the normal coordinates is just the product of 
two one-dimensional harmònic oscillators, whereas the formulation in the original coordinates 
is somewhat more complicated, and features the interdependence of the two oscillations of the 
respective single pendula. Limiting cases, e.g. m\ 3> rri2, are easily recovered yielding two 
(almost) independent oscillations (r?^ 2 'g<fi^>2 becomes negligible). For the sake of simplicity I 
have considered only the case l\ = I2 = l- A generalization to l\ 7^ I2 is easily done, and left to the 
reader. The generalization to higher dimensions is also straightforward; again we would obtain 
a product of n-dimensional harmònic oscillators in terms of normal coordinates, i.e., £ = Cq, 
respectively q = C -1 £, with C the corresponding affine transformation. 

These two examples thus show in a nice way the possibility of solving path integrals for 
kinetic energy couplings. Other interesting cases would be a particle in rotating magnètic fields, 
an arrangement important for trapping electrons, and non-orthogonal coordinate systems in 
spaces of constant curvature [18|. This will be studied in subsequent investigations. 
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